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Abstract 

In this paper, we improve the calculation of the relic gravitational waves (RGW) in two aspects: First, we 
investigate the transfer function after considering the redshift-suppression effect, the accelerating expansion 
effect, the damping effect of free-streaming relativistic particles, and the damping effect of cosmic phase 
transition, and give a simple approximate analytic expression, which clearly illustrates the dependent rela- 
tions on the cosmological parameters. Second, we develop a numerical method to calculate the primordial 
power spectrum of RGW at a very wide frequency range, where the observed constraints on n s (the scalar 
spectral index) and Ps(ko) (the amplitude of primordial scalar spectrum) and the Hamilton- Jacobi equation 
are used. This method is applied to two kinds of inflationary models, which all satisfy the current constraints 
on n s , a (the running of n s ) and r (the tensor-scalar ratio). We plot them in the r — £l g diagram, where Q g 
is the strength of RGW, and study their detection by the CMB experiments and laser interferometers. 
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I. INTRODUCTION 

Recently, a lot of observations on the CMB power spectra [l-^] and the large scale structure 
(LSS) [4] have supported inflation as the good phenomenological model to describe the evolution of 
the universe at very early stage, which naturally answers the origin of the primordial fluctuations 
with a nearly scale-invariant and gaussian spectrum. In addition to the density perturbations, 
inflationary models also predict a stochastic background of relic gravitational waves (RGW), which 
is also called the tensor perturbations. The amplitude of RGW directly relates to the energy scale 
of inflation. Although this background has not yet been observed, and only some loose constraints 
have been achieved [3, Q], but its detection would provide incontrovertible evidence that inflation 
actually occurred and would set strong constraints on the dynamic of inflation [6] . So it is always 
regarded as the "smoking-gun" evidence for inflation. 

There are mainly two kinds of experiments are underway to detect the RGW at different 
frequency. One is the CMB experiments, which can find the RGW by observing the CMB B- 
polarization power spectrum [7]. This method is sensitive to the waves with very low frequency, 
v e (KT 17 , lCT 15 )Hz. Now, the first-three -year results of WMAP [2] have not found the evidence 
of the gravitational waves, and only give a constraint r < 0.28(95% C.L.), where r is the so-called 
tensor-scalar ratio. The next experiment, the Planck satellite [Sj, has higher sensitivity to polar- 
ization, which is scheduled for launch in 2007, and expects to observe the RGW if r > 0.1. The 

Q 

ground-based experiment, Clover (Cl-Observer) is also under development [9j, which is expected 
to observed the RGW if r > 0.005. Another kind of important experiments are the laser interfer- 
ometers, including the BBO (Big Bang Observer) and DECIGO (DECihertz Interferometer 
Gravitational wave Observatory) 111], which can detect the gravitational waves with very high 



frequency v ~ 0.1Hz. The former can detect the RGW when Q, g > 2 x 10~ 17 is held, where Q 9 
is the strength of RGW at 0.1Hz, and the latter expects to observe if £2 g h 2 > 10~ 20 . It should 
notice that the waves with very high frequencies can be observed by the electromagnetic resonant 
systemfljj]. This is also an important method to detect the relic gravitational waves. 



A lot of works have been done to study the RGW detection by these experiments [13j|. In the 
previous work we have discussed the predicted values of RGW (r and Q g ) for some kinds of 
inflationary models, where we have used a simple power-law function to describe the primordial 
power spectrum of RGW, which is a very good approximation for the waves with very low frequency, 
but for the waves with high frequency, this may generate large error. At that work, we have not 
considered the damping effect of cosmic phase transition on the RGW, such as the QCD transition 
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15l4l7l|. e + e _ annihilation and so on. In this paper, we will discuss this topic more precisely: 
First, we will consider the damping effect of a general cosmic phase transformation, which can 
been described by a simple damping factor. And then we will give a simple form of the total 
transfer function, which applies to the waves with v 10 _16 Hz. This function is dependent on the 
values of Q m , the present energy densities of vacuum and matter, respectively; the value of to, 
the age of the universe; the value of Hq, the present Hubble constant; the values of g* and g* s , the 
effective number of relativistic degrees of freedom when the waves exactly crossed the horizon; and 
the fraction / of the background (critical) energy density of the free-streaming relativistic particles 
in the universe when the waves exactly crossed the horizon. So this function includes abundant 
cosmic information. Second, we will use a numerical method to calculate the primordial power 
spectrum of RGW, where the Hamilton-Jacobi formula is used. Compared with the simple power- 
law form, this numerical result has little change for the value of £l g when r is smaller, r < 0.02. 
But when the value of r is larger, the numerical result is obviously smaller than that of the simple 
power-law approximation. 

The organization of this paper is as follows: in section 2, we will simply review the evolutive 
equation of the RGW. In section 3, we mainly discuss the damping effects. In section 4, we will 
introduce the numerical method by discussing two kinds of inflationary models. At last, we will 
give a conclusion and discussion in section 5. 



II. THE RELIC GRAVITATIONAL WAVES 



Incorporating the perturbation to the spatially flat Robertson- Walker (FRW) spacetime, the 
metric is 

ds 2 = a(T) 2 [dT 2 - (Si, + hifi&jdx?] , (1) 

where a is the scale factor of the universe, r is the conformal time, which relates to the cosmic time 
by adr = dt. The perturbation of spacetime hij is a 3 x 3 symmetric matrix. The gravitational 
wave field is the tensorial portion of hij, which is transverse-traceless dih 13 = 0, 5 lJ hij = 0. Since 
the gravitational waves are very weak, \hij\ <C 1, one needs just study the linearized evolutive 
equation: 

d^V^d^hij) = 167rGa 2 (T)ny , (2) 

where Iljj is the tensor part of the anisotropy stress, which satisfies Tin = and djlij = 0, and 
couples to h-ij like an external source in this equation, which is always generated by the free- 
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20l | . and so on. It is convenient to 



streaming relativistic particles [la . [19(, the cosmic magnetic 
Fourier transform as follows 

Mr,x) = £ v^G | ^^4 A )(k)^(r)e ikx , (3) 

A 

n,,(r,x) = £ VlfoFG I ^^4 A) (k)n^(r)e ikx , (4) 

A 

where A = ' + ' or "x" labels the two polarization states of the gravitational waves. The polarization 
tensors are symmetry, transverse-traceless fe*e^ (k) = 0, 5* J e^ (k) = 0, and satisfy the conditions 
e ( A )*j(k)e^ ^(k) = 25\\i and e^(— k) = (k). Since the RGW we will consider is isotropy, and 
each polarization state is same, we can denote h^\r) by %(r), and EE^ (r) by Ilfc(r), where 
k = |k| is the wavenumber of the gravitational waves, which relates to the frequency by v = k/2n 
(the present scale factor is set ao = 1). So Eq.([2]) can be rewritten as 

h k + 2-h k + k 2 h k = 16vrGa 2 (r)n fe (r) , (5) 

where the overdot indicates a conformal time derivative d/dr. 

The RGW generated during the early inflation stage. Inflation is an extremely attractive idea 
to describe the very early universe, which has received strong support from the observations of 
CMB anisotropies and from studies of the large-scale distribution of galaxy. In this paper, we will 
only consider the simplest single field models. This kind of models is enough to account for the 
current observations on n s , a, and r. In the context of slow-roll inflationary models, the most 
observables depend on three slow-roll parameters 21] 
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where Mp\ = (87rG) -1 / 2 = mp\/\^87r is the reduced Planck energy. In the following discussion, I 
will use the unit Mp\ = 1 and mp\ = y/8ir. V((j)) is the inflationary potential, and prime denotes 
derivatives with respect to the field (j). Here, ty quantifies "steepness" of the slope of the potential, 
r]v quantifies "curvature" of the potential and £y quantifies the "jerk". All parameters must be 
smaller than one for inflation to occur. The most important prediction of the inflationary models 
is the primordial scalar perturbation power spectrum, which is nearly gaussian and nearly scale- 
invariant. This spectrum is always written in the form 

/ l. \ n,(*o)-l+|ala(*/*o) 

Ps(k) = P s (k ) -) " , (7) 
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where a = dn s /dlnk, and ko is some pivot wavenumber. In this paper, k = 0.002MPC" 1 is used. 
The observations of WMAP find P s (k ) ~ 2.95 x 1(T 9 A(k ) and A(k ) = 0.813l{j;g^ [2]. Another 
key prediction of inflationary models is that the existence of the RGW. The primordial power 
spectrum of RGW is defined by 

P T (k) EE ~^Kh k . (8) 

The strength of the gravitational waves is characterized by the gravitational waves energy spectrum 



where p c = 3Hq/8ttG is the critical density and Hq = WOh km s 1 Mpc 1 is the 



constant. One can relate Q g to the primordial power spectrum by the formula 

1 



present Hubble 



19] 



fy(*0 = J^2k 2 P T (k)Tf(k) , (10) 



where the transfer function TAk) reflects the damping effect of the gravitational waves when 
evolving in the expansion universe. It is convenient to defined a function T(k) = k 2 T 2 /12Hq, so 
the strength of RGW becomes £l g (k) = T(k)Pr(k). In the following sections, we will discuss T(k) 
and Pr(k), respectively. 



III. THE TRANSFER FUNCTION 



In this section, we will discuss three kinds of damping effects: First we ignore the anisotropy 
stress in Eq.([5]), and only consider the redshift-suppression effect. So Eq.([5]) becomes 

h k + 2-h k + k 2 h k = . (11) 
a 

This is the evolutive equation of RGW in vacuum, which only depends on the evolution of the 
scale factor a(r). It is clear that, the mode function of the gravitational waves behaves simply in 
two regimes when evolving in the universe: far outside the horizon (k <C aH), and far inside the 
horizon (k 3> aH). When waves are far outside the horizon, the amplitude of h k keeps constant, 
and when inside the horizon, they damp with the expansion of the universe 

. (12) 

In the simple cosmic model, the evolution of the universe can be separated into three stages: the 
radiation-dominant stage, the matter-dominant stage, and the vacuum-dominant stage. In this 
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model, by numerically integrating the Eq. ljlip . one found the transfer function can be approxi- 
mately described with a damping function (for the waves with k 3> 10~ 18 Hz) Q, H . 23]: 



™ = vk>mr +1 *(£) + ™(£) • (13) 

where k eq = 0.073fi m /i 2 Mpc _1 is the wavenumber corresponding to the Hubble radius at the time 
that matter and radiation have equal energy densities. And tq = 1.41 x 10 4 Mpc is the present 
conformal time. fi m and £l\ are the present energy densities of matter and vacuum, respectively. 
It is obvious that, when k <C k eq , which entered the horizon in the matter-dominant or vacuum- 
dominant stage, t\{k) ex k~ 2 , but when k ^> k eq , which entered the horizon in the radiation- 
dominant stage, ti(k) oc k^ 1 , which is for the different evolution of scale factor in different stages. 
The factor Q JV J^l\ is the effect of accelerating expansion, which has been discussed in the previous 



works 



23, 



2J. 



The second is the damping effect of the free-streaming relativistic particles [18(] , especially the 
neutrino, which can generate the anisotropic stress 11/% on the right-hand of the Eq.©, when it is 
the free-streaming relativistic particles. This effect was first considered by Weinberg, where the 
Eq.© can be rewritten as a fairly simple integro-differential equation. The solution shows that 
anisotropy stress can reduce the amplitude for the wavelengths that re-enter the horizon during 
the radiation-dominated stage, and the damping factor is only dependent on the fraction / of the 
background (critical) energy density of the free-streaming relativistic particles in the universe. The 
effect is less for the wavelengths that enter the horizon at later time. A lot of works have been 
done to simplify this effect, and in Ref.[3], the authors found it can be approximately described 
by a transfer function t2 for the waves with v > 10 _16 Hz (which re-enter the horizon at the 
radiation-dominant stage), 

_ 15(14406/ 4 - 55770/ 3 + 3152975/ 2 - 48118000/ + 324135000) 
h ~ 343(15 + 4/)(50 + 4/)(105 + 4/)(180 + 4/) ' ( ' 

When the wave modes (10 _16 Hz< v < 10~ 10 Hz) re-enter the horizon, the temperature in the 
universe is relatively low (< lMeV), we are fairly confident that the neutrino is the only free- 
streaming relativistic particle. So we choose / = 0.4052, corresponding to 3 standard neutrino 
species, the damping factor is 0.80313. But for the waves with very high frequency (v > 10~ 10 Hz), 
the temperature of the universe is very high when they re-enter the horizon, and the value of / is 
much uncertain. Thus, the detection of RGW at this frequency offers the probability of learning 
about the free-streaming fraction / in the very early universe. 
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The third is the effect due to the successive changes in the relativistic degrees of freedom during 
the radiation-dominant stage, here we also call it the effect due to the cosmic phase transition, 
which includes the QCD transition, the e + e~ annihilation, the electroweak phase transition and 
so on. In an adiabatic system, the entropy per unit comoving volume must be conserved 
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25|, 



27T 2 

S(T) = s{T)a 3 (T) = constant, and s(T) = g* s (T)T 3 , (15) 

45 

where the entropy density, s(T), is given by the energy density and pressure; s = (p + p)/T. Com- 
bining it with the expressions of energy density and pressure in the radiation-dominant universe, 

p{T) = ^(T)T\ P (T)= l -p(T). (16) 

one can immediately get the relation 

-4/3 _4 / iv x 

p oc g*g* s a , (17) 

where we have defined the "effective number of relativistic degrees of freedom", and g* s , following 



the Ref. 



15|, 



251 ] . These quantities, g* and g*s, count the effective number of relativistic species 
contributing to the radiation energy density and entropy, respectively. From this relation, one can 
find that, if the phase transitions are not considered, g* and g* s are all constant, and this relation 
becomes to the general case of p oc a -4 . However, it does not always hold, as some particles would 
become non-relativistic before the others and stop contributing to the radiation energy density. 
In other words, the evolution of p during the radiation era is sensitive to how many relativistic 
species the universe has at a given epoch. As the wave equation of gravitational waves constraints 
(a/a)hk, the solution of hk would be affected by and g*s via the Friedmann equation: 

H(r) \ 2 _{g*\ fg* s \ ~ 4/3 Qr (±y 4 + Qm ( ±) " 3 + o Aj (is) 

ic 



-no / \g*oj \g*so/ \ a o/ \ a o 

where the subscript denotes the quantities with the present values. Here we have considered t 



1. 



Friedmann equation in the ACDM universe, which is supported by a number of observations 
0, 26]. Inserting this into Eq. ljlip . one can numerically calculate the value of hk [3], which 
would spend a very long computer time, since one must integrate that equation from the end of 
the inflation to present time, and calculate the waves from v = 10~ 16 Hz to 0.1Hz which we are 
interested. Here we will give an approximate method, which can describe this effect by a simple 
factor £3. We consider the wave hk with the wavenumber k, which crossed the horizon at a = a^, 
and the corresponding Hubble parameter is Hk- So one has k = akHk/ao. One knows that,when 
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the waves are in the horizon, /ifc(r) oc l/o(r), damping with the expansion of universe, and when 
the waves is out the horizon, the = constant, keeping its initial value. So one can define a factor 

F k = ^ = ^ . (19) 

where is the conformal time at the beginning of the radiation era. During the radiation era, one 
has 

/ \ 1/2 / \ -2/3 / \ -2 

H = B (<t-) (*=-) (±) , (20) 

\g*oJ \g*soJ \aoj 

\i<2 

where B = HqQ, t , is a constant. Using the Eqs. (|19p . (|20p and the relation of k = df-H^/ao, one 
gets 

Ft = ^(9jm) 1,2 (s^MY 2, \ (21) 

k V 9*0 ) V 9*s0 ) 
where T)% is the temperature when the wave exactly crossed the horizon. First we can assume 
g* = 5*o an d g* s = g*s0 are always satisfied, which is the condition without changes in the rela- 
tivistic degrees of freedom during the radiation era, which follows that = B/k. Inserting this 
into Eq. fllOj) . one finds the value of £l g (k) is independent of the wavenumber k. However, here 
we are interested in the condition with changing and g* s , and the factor £3 exactly denotes the 
difference of these two conditions, i.e. 

t3 = ^ = (^y 2 ( s _^itiy 2, \ (22) 



Fk \ g*o J V g*so 

where 5*0 = 3.3626 and g* S Q = 3.9091. This factor depends on the values of <?* and g*s at the early 
universe. Fig.[l], presents the evolution of the values of and g* s , which shows that the value 
of has an obvious accretion when T > O.lMeV. The difference between g* and g* s only exists 
when T < O.lMeV. In the expression of Q g (k), this effect is described by a factor i 3 2 . Compared 
with the accurately numerical calculation, this approximation has the error smaller that 10%. The 
total transfer function is the combination of these three effects 

Tf(k) =tixt 2 xt 3 , (23) 

where t\ is most important, which approximately shows the evolution of RGW in the expanding 
universe. The function of ti has the most uncertain in this discussion. In the extreme condition 
with / = 0, t2 = 1 is held, i.e. no damping; and in another extreme condition with / = 1, t 2 = 0.59 
is held, this function arrives at its smallest value. In the case of / = 0.4052, t% = 0.80321 only 
contributes a damping factor 0.645 for the strength of the RGW. The value of £3 is fairly small. 



For the extreme condition with T k > 10 6 MeV (k > 2 x 10 _4 Hz), one has g* — g* s — 106.75 in the 
Standard Model (3* = g* s = 228.75 in the MSSM), t 3 = 0.62 (t 3 = 0.55 in MSSM) only contributes 
a damping factor 0.38 (0.30 in MSSM) to the strength of the RGW. 

The experiments which can directly detect the RGW are all sensitive to the waves with k 3> k eq , 
which re-entered the horizon during the radiation era. From the previous discussion, one can get 
a simple expression of all these damping effects: 

j-ys _ ( 15 \ (VL m \ 2 ( g*(T k )\ ( g*s(Tk)\~ 4/S 



x 



SklqH^ J \ Q A J V 9*0 J \ g*s0 



( 15(14406/£ - 55770/| + 3152975 /| - 48118000/ fc + 324135000) \ 2 

V 343(15 + 4/ fc )(50 + 4/ fc )(105 + 4/ fc )(180 + 4f k ) ) " ( } 



where f k is the value of the function / when wave h k exactly crossed the horizon. This function is 
dependent on the values of CI a, the present energy densities of vacuum and matter, respectively; 
the value of To, the age of the universe; the value of Hq, the present Hubble constant; the values of 
g*{T k ) and g* s (T k ), the effective number of relativistic degrees of freedom; and the fraction f k of 
the background (critical) energy density of the free-streaming relativistic particles in the universe. 
So this function includes abundant cosmic information. Using this, the strength of RGW becomes 

n g (k) = P T (k)T(k) . (25) 

Here we are interested in the waves with u = 0.1Hz, which is the sensitive frequency of laser 
interferometers, BBO and DECIGO. Choosing the cosmic parameters h = 0.72, Q m = 0.27, Ua = 
0.73, g* = g* s = 106.75 and f k = 0, one gets 

T(k) = 4.15 x 10~ 7 , Q g (k) = 4.15 x 10~ 7 P T (A;) . (26) 

IV. THE PRIMORDIAL POWER SPECTRUM OF RGW 

The primordial spectrum of RGW is always described in a simple form 

(, \ nt(ka)+\a t \Ti(k/ko) 
¥ J ' , (27) 

where n t (k) is the tensor spectral index, and at = dn t /d\nk is its running. In the single- field 
inflationary models, a standard slow-roll analysis gives the below relations of observable quantities 
and slow-roll parameters, 



r 

n t = -77 > a t 



r 



8 16 

(n s -l) + - , r= -(l-n s ) + — W , (28) 



where r(k) = Pr(k)/Ps(k), is the so-called tensor-scalar ratio. These formulae relate rit and at 
to the other two functions n s and r, which are easily to observe. But the relation between r and 
n s is dependent on rjy, which depends on the specific inflationary potential. Inserting these into 
Eq.([27|), one gets 



k 



-i+^[("»-i)+S]M*/*o) 



P T (k) = Ps(k ) x r x , (29) 

where r denotes the tensor-scalar ratio at k = ko, i.e. r = r(ko) , which is also held in the following 
sections. So the primordial spectrum of RGW only depends on n s and r. The recent constraints 
come from the observations of three- year WMAP 2;] , which are 

n s = 0.951+^ (68% C.L.) , r < 0.28 (95% C.L.) . (30) 

Using the Eq.(f2"B"]). one gets 



/ k \ -i+w[( n «- 1 )+S] ln ( fc /*o) 

Q g (k) = 9.98 x HT 16 r ( — J , (31) 

where we have chosen A{ko) = 0.813. We have plotted the function Vt g (denoting Q g = O g (/ei), 
and k\ = 0.1Hz) dependent on r in Fig. [2], where n s = 0.951 is used. This result is consistent with 
our previous work, a larger r leads to a larger Q, g . We well know that the formula in Eq. (|27[) is a 
very good approximation when the wavenumber k is not much larger (or smaller) than fco- But it 
may be not a good approximation at k%, which is more than 16 order that &o- So it is necessary 
to numerically calculate Pr(k), but it is not easy for the exactly numerically calculation, since one 
must calculate the spectrum in a very wide range in wavenumber (larger than 16 order), and for 
each k, one must integrate it from the initial condition to the end of the inflation. In this section, 
we will use a semi-numerical method to calculate the primordial power spectrum of RGW. We will 
introduce this method by discussing two kinds of inflationary models, which all satisfy the current 
constraints of n s , a and r. 

First we consider the model with potential (Mod. 1.1) V(<p) = A 4 (^///) 2 , which belongs to the 
large- field model, and predicts a fairly larger r 2?|]. From Eq.([6]), one gets 

2 2 

where we have denoted Mp\ = 1. So the slow-roll condition requires that <p 3> v2 ; the so-called 
large-field model. At the end of inflation, ey = 1 is satisfied, which leads to (p e nd = v2- In the 
initial condition, one has 21] 



V 

n s -l = -6e v + 2r) V , r = 16ey, Ps(/co) = j^l^J ' ( 33 ) 
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which follows that 

4>ini= V8/(l~n s ), r = 4(l-n s ) A 4 //x 2 = 0.75ir 2 P s (k )(l - n s f . (34) 
Inserting these into the Hamilton-Jacobi formula, 

2[H'{4>)] 2 -?,H\ ( p) = -V{<t>) , (35) 



one can immediately get the function H{(p) by the numerical calculation. We define the e-folds 
a = ciini = 1 i.e. N = 0. The relation of N and <j) is 



number N, and the scale factor is a = a,j„ie N ■ When kn crossed the horizon, we set the scale factor 



dN = - 2 Jl> (36) 

where H is the Hubble parameter during inflation, and H' = dH/d(f). One can define a Hubble 
slow-roll parameter e = 2{H' / H) 2 , so the primordial power spectrum of RGW is (to the first 



slow-roll order) 



Pj (7,'j = - 2 



H 2 



(37) 

k=aH 



where c = 4(ln 2 + 7) — 5 ~ 0.0814514 (with 7 the Euler-Mascheroni constant) is a constant. Using 
the Eqs. (j35|) and (p6|) . one can numerically calculate H(N). Inserting it into Eq. (p7|) . one can get 
the primordial spectrum of RGW, at the same time the total e-folds iV is also got. The value of 0, g 
is also got by using the Eq. (|26|) . We have plotted the values of r and £l g in Fig. [2], where we have 
chosen n s E [0.94, 0.98]. It is easily found the value of r is in the range r £ [0.08, 0.24]. Compared 
with the approximate formula (|3ip . the numerical value is much smaller: When n s = 0.951, the 
value is only one third of the approximate value. 

Then let us consider another model V(<p) = A 4 [l — (cj)/fi) 2 ], which belongs to the small-field 
model, and predicts a very small r (2^]. From Eq.©, one gets 

1 [ 2x1 jx I 2 2/y 2 

where x = (J)/[i. At the end of inflation, 4> en d = ^ i.e. x enc [ = 1, where V = is satisfied. The 
initial value of x must be very small to account for the slow-roll condition. Since it can not be got 
from the observed n s and Ps(ko), we must set it before the calculation. First let us consider the 
model with x ini = 0.1 (Mod. 2.1), using the Eq, ([3"3"j) . one immediately gets 

H 2 = 4.89746/(1 -n s ), r = 0.06667(1 - n 8 ), A 4 = 0.99693(1 - n s )P s (k ); (39) 
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Second W6 consider th.6 model with x% n i — 0.2 (Mod. 2. 2), which follows that 

^ = 6.07693/(1 -n s ), r = 0.22857(1 - n s ), A 4 = 3.52486(1 - n s )Ps(ko); (40) 

The third model has Xj n j = 0.3 (Mod. 2. 3), which follows that 

l? = 7.72854/(1 - n 8 ), r = 0.45(1 - n s ), A 4 = 7.32086(1 - n s )P s (k ). (41) 

Then using the Hamilton- Jacobi formula in Eq. (j35p and the relation of N and (ft, here it becoming 
d(j)JdN = 2H' /H, one can also get the function H(N). Using the formula (|37p . the values of r, 
Pr(k), £l g (k) and N are also got, which have been plot in Figs. [2] and [3]. From Fig. [3], one finds a 
larger n s leads to a larger N, which is held for all these four inflationary models. When n s = 0.951, 
N = 41.96 for the Mod.1.1, and N = 62.47 for the Mod. 2. 3, which are in the region of N G [40, 70], 
and acceptable But for the Mod.2.1, N = 97.90, and for the Mod. 2. 2, N = 74.59, which are 
too large to suitable. From Fig. [2], one finds that when n s 6 [0.94, 0.98], r < 0.02 is satisfied for 
the Mod.2.1, Mod. 2. 2, Mod. 2. 3, the very small values. And values of £l g are exactly same with the 
approximate results. So one gets a conclusion: when r is small, the formula (|3ip is a very good 
approximation, but when r is larger (r > 0.1), the approximate formula (|3ip is not very good, and 
the numerical calculation is necessary. 

V. CONCLUSION AND DISCUSSION 

Inflation has received strong supports from the observations of the CMB and LSS. As a key 
prediction of inflationary models, the detection of RGW would provide incontrovertible evidence 
that inflation actually occurred and set strong constraints on the dynamics of inflation. A lot 
of experiments are under development for the RGW detection, which mainly include two kinds: 
The CMB experiments, including Planck, Clover, and others, and the laser interferometers, in- 
cluding BBO, DECIGO and so on. For investigating the detection abilities of these two kinds of 
experiments, it is convenient to study the distribution of the inflationary models in the r — £l g 
diagram. So it is necessary to accurately calculate the RGW at all frequency range. In this paper, 
we improved the previous calculation in two aspects: First, we studied the transfer function after 
considering the redshift-suppression effect, the accelerating expansion effect, the damping effect of 
free-streaming relativistic particles, and the damping effect of cosmic phase transition, and gave 
a simple approximate formula of the transfer function, which applies to the waves with k > k eq . 
This function depends on the values of the cosmic parameters: TO , Hq, k eq , To, g*, g* s , and 
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fk- Second, we have developed a numerical method to calculate the primordial power spectrum 
of RGW, especially at high frequency, where the observed constraints on n s and Ps(ko) and the 
Hamilton-Jacobi equation are used. We applied this method to two kinds of inflationary models, 
which all satisfy the current constraints on n s , a and r. 

From Fig. [3], one can find: in all these inflationary models, a larger n s follows a larger N. For 
the first kind of model, when n s > 0.97, the value of N > 70 is satisfied, which is unsuitable. To 
account for the constraint of N £ [40, 70] , n s can only be in a very narrow region n s £ [0.948, 0.970] . 
For the second kind of model, the initial conditions of Xi n i = 0.1 and 

x ini — 0.2 are not acceptable, 

which predict too large e-folds. The condition of Xi n i = 0.3 is suitable, which predicts N = 62.47 
when n s = 0.951. But to account for the constraint N < 70, n s < 0.956 must be satisfied. From 
Fig. [2], one found that, for the Mod. 1.1, when n s £ [0.94, 0.98], the value of r is in the region 
r £ [0.08, 0.24], which are mostly in the sensitive region of Planck satellite. The value of is in 
the region of [5.6xl0~ 17 , 2.2 xl0~ 18 ]. In the most region of n s , a larger n s follows a smaller r, and 
corresponds to a larger Q g , which is an unexpected result. This is an obvious difference from the 
result of the approximate formula. When n s = 0.951, Q g = 1.3 x 10" 17 , which is in the sensitive 
region of ultimate DECIGO, but beyond the sensitive limit of BBO. This value is only one third 
of the value from the approximate formula in Eq. (|3ip . For the Mod. 2.1, Mod. 2. 2 and Mod. 2. 3, 
a larger n s follows a smaller r and a smaller Q g . These models predicted a very small r, when 
n s £ [0.94, 0.98], r < 0.02 is always satisfied, and the value of £l g is exactly same with the the 
value from the approximate formula in Eq. (|3ip , For the Mod. 2. 3, which predicted an acceptable 
e-folds, the values of r are all in the sensitive region of Clover, but beyond which of Planck; the 
values of £l g are all in the sensitive region of ultimate DECIGO, but beyond which of BBO. 
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